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In this invited presentation I review some recent developments in the theory of charmonium that 
appear likely to be of importance for future experimental studies in this field. The specific areas con- 
sidered are double charmonium production, LQCD studies of charmonium, recent results for hadron 
loops, cc production cross sections at PANDA, charm molecules, and two recent developments, 
"charmiscelleny" . 



Introduction 

' Charmed hadron spectroscopy has undergone a renaissance in recent years, due to an impressive increase in exper- 
imental results, which are often in striking disagreement with theoretical expectations regarding these hadrons. This 
$H renaissance began with the discovery of the surprisingly light charm-strange mesons [1, 2], and has continued with the 
observation of the charm-meson molecule candidate X(3872) and a "zoo" of latter-alphabet X,Y and Z states, which 
include conventional 2P cc candidates, hybrid charmonium candidates, and "charged charmonia" , which are of course 
charmed meson molecule candidates. For recent reviews of the experimental status of this field, see Refs.fs', 4^. Many 
of these exciting experimental discoveries have rather surprisingly originated from facilities that were not originally 
intended for research in this area, notably the e"'"e~ "B factories" , which have proven to be copious sources of charmed 
hadrons. 

■ In this report I will discuss recent developments in our understanding of the theory that underlies charmonium 
states and related topics. My criteria in selecting material for this talk are that the results should be clearly relevant 
i' to experiment, reasonably recent, and likely to be relevant to future experimental studies in this field. The specific 
O , topics chosen for review are as follows: 
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1. Double charmonium production 

2. LQCD studies of charmonium 

3. Recent results for hadron loops 

4. Charmonium production cross sections at PANDA 

5. Charm molecules 
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' 6. "Charmiscelleny" ; two recent developments 



I. DOUBLE CHARMONIUM PRODUCTION 

One of the interesting new experimental techniques that has been exploited recently in charmonium studies is 
"double charmonium production", in which e''"e~ annihilation produces two charmonium resonances (see Fig[T]). In 
the results reported thus far, a J/^p has been used as a trigger; the missing mass of the recoiling system shows clear 
peaks at the masses of the rjc, xo a-nd rj'^, and a higher- mass enhancement referred to as the X(3943) is also observed. 
(See for example results reported by the Belle (^Q and BABAR 7] collaborations.) 

The double charmonium production process is especially interesting in that it gives us access to J^^ ^ 1 cc 
states at e"'"e~ machines. This includes C = (-I-) cc mesons in particular, which are otherwise rather difficult to 
produce; typically they have been studied using hadronic production, as in pp annihilation at Fermilab, in radiative 
transitions from higher-lying J^*^ = 1 states, or in the two-photon process e^e^ — > 6+6^77,77 cc, which has 
rather small cross sections. 
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Although Figl2] shows that this is clearly a successful approach for producing C = (+) charmonia in e+e" anni- 
hilation, it would be even more interesting if the double charmonium production mechanism were well understood 
theoretically; in this case one could ideally predict the production cross sections for various cc states, and non-cc can- 
didates could be identified through their anomalous cross sections. In addition one might hope to use this mechanism 
to observe the as yet unreported — 2^ cc D-wave states ric2 and hc2', these are expected to be quite narrow, due 
to the absence of open-charm hadronic decay modes. 




C={+) cc 



FIG. 1: Double charmonium production in e^e annihilation. 
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FIG. 2: States seen in double charmonium production (recoiling against the J/V') by Belle 



The theoretical mechanism of double charmonium production has been studied in some detail by three groups, 
at ANL, Beijing and Oxford. In their important initial study the ANL group (Braaten and Lee Q) assumed a 
nonrelativistic one-gluon-exchange (OGE) mechanism for production of the "second" cc pair. (One might a priori have 
anticipated important contributions either from pQCD-type mechanisms such as OGE, or from the nonperturbative 
decay mechanisms that appear to dominate light hadron strong decays, as are described by the ^Pq model.) Braaten 
and Lee concluded that the OGE decay mechanism gave a characteristic large amplitude for e+e" J/ip + Xco 
relative to the other Xcj states, with approximate relative cross sections of 

a(e+e- ^ J/V' + (Xc2 : Xci : Xco)) « 3 : 2 : 10 . (1) 

This explains the dominance of the XcO over the other IP states that is observed in the data (see for example Figl2|). 
However it was also noted that these leading order (OGE) NR pQCD cross sections are about an order of magnitude 
smaller than the experimental results. 

As this was a rather complicated pQCD calculation, it was clearly important to have an independent check of these 
results. This work was confirmed and extended in a series of papers by the Beijing group of Chao, Liu, He, Ma and 
Zhang; see Ref.Q and references cited therein. Liu, He and Chao (oj have recently reported very interesting predictions 
of various double charmonium cross sections, including results for states that have not yet been observed, such as 
the narrow ?7c2 and hc2 states. Their results for these cross sections (relative to the reference channel J/V' + Vc) 
are shown in Table HI evidently these predictions imply that the double charmonium production mechanism has 
considerable potential for the discovery of new states, given an improvement over current statistics by one to two 
orders of magnitude. 

Regarding the problematic overall scale of these cross sections, Zhang, Ma and Chao Ji| considered NLO corrections 
to the OGE model, and found that these considerably reduced the reported discrepancy between the observed and 
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TABLE I: Double charmonium production cross sections relative to J/tp + rjc at ^/s = 10.6 GeV found by Liu et al. 

calculated double charmonium production cross sections. Thus the well-known problem of the scale of these cross 
sections may have been due to the neglect of higher-order pQCD effects. 

The Oxford group of Close and Downum [11] has also investigated double charmonium production. They too 
confirmed the NR pQCD OGE calculation of Braaten and Lee analytically, and have in addition considered nonper- 
turbative decay mechanisms of ^Pq type (see Burns, Close and Thomas [12,]). Close and Downum found that in the 
latter case the final state J/ipXcO would no longer be the dominant J/tpXcj niode, which argues in favor of pQCD 
processes for double charmonium production. They note that in contrast to the double charmonium results, in the 
light meson sector one sees vector -|- tensor dominance, as in e"'"e~ — )■ 10/2, which is a prediction of the nonperturbative 
decay mechanism. 

II. LQCD CHARMONIUM 

Here I will be very succinct, as the very important new work in this area will be discussed in detail in a dedicated 
plenary talk by C.Thomas [13]. 

Although LQCD studies of charmonium spectroscopy have existed in the literature for some time (see for example 
Ref.[l3]), these studies have not had especially high impact because they have generally confirmed the predictions of 
cc potential models and experimental results for the known spectrum of cc states. More dramatic early predictions 
from LQCD have included mass estimates for quarkonium hybrids, such as ca. 4.4 GeV for the lightest 1 ^ exotic 
charmonium hybrid. 

Recently however very impressive new results on charmonium from LQCD have been rep orted, notably by the 
JLAB lattice QCD collaboration of Dudek, Edwards, Peardon, Richards and Thomas [1^, [l^. The most interesting 
results for charmonium involve determinations of the photocouplings of charmonium sector resonances, including 
J^'^-exotic charmonium hybrids jl5| . The importance of these new results for the JLAB experimental program would 
be hard to overstate. The 12 GeV upgrade at JLAB and the Hall-D GlueX experiment in particular are predicated 
on the hope that the photocouplings of hybrid mesons are sufficiently large for the spectrum of these states to be 
clearly established through a high-statistics study of meson photoproduction. Although models of hybrids suggest 
that their photocouplings should not be small, these of course make assumptions regarding the nature of hybrids that 
are conjectural, and may be inaccurate. Without strong theoretical evidence of large hybrid photocouplings, especially 
for the "smoking-gun" J^'^-exotics, the plan to establish these states experimentally using photoproduction at JLAB 
is not well justified theoretically. Clearly, lattice QCD studies of the photocouplings of hybrids have an extremely 
important role to play in establishing the sensitivities required for a definitive study of hybrids using photoproduction 
reactions, due to their better controlled systematic uncertainties. 
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TABLE IL Ml radiative decay widths for charmonium and related states from LQCD, as reported by Dudek et al. [Tsl ]. 
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An example of the very important new results reported by the JLAB LQCD group, from the viewpoint of using 
photoproduction to identify hybrids, appears in their table of predictions for Ml radiative partial widths (Table lll| . 
First note the reasonably accurate result for the "ground state" (15 — )• IS) Ml transition J/ip — >• rjcj; lattice data and 
experiment are shown together in the rightmost column, and Tthy./^expt. = 1.36±0.22. This shows that Ml transitions 
are reasonably well reproduced in this LQCD study. The less significant results for 25 — > 15 V'' — >■ Vcl involve a well 
known difficulty; since the naive quark model matrix element is zero due to orthogonal spatial wavefunctions, this is an 
attempt to extract an intrinsically small amplitude. The most interesting result here is the corresponding Ml partial 
width predicted for the non-exotic J^'~^ = 1~~ charmonium hybrid ''Yhyb.''\ which is r{Yhyb. — '7c7) = 42 ± 18 keV. 
By the standards of Ml transitions within cc states this is a large rate, and implies that radiative transitions between 
hybrids and ordinary quarkonium ground states are not significantly suppressed. Even more interesting is the Ml 
radiative decay width reported by this reference for the "classic" hybrid exotic with J^'^ — 1 ^, which is 

r(r;ci(l^+ exotic) ^ J/V'7) - 100 keV. (2) 

This is a very large partial width for an Ml transition, and is comparable to the known El partial widths for transitions 
between members of low-lying cc multiplets, such as iJj' — ^ XcJ7 a-nd Xc.i — ^ J li^l- This is exciting indeed, as it implies 
that the GlueX program to search for exotic hybrids using photoproduction at JLAB is well motivated. This result 
may also imply that largely hybrid states (at least with heavier quarks) can be identified through their radiative 
partial widths, such as anomalously large Ml widths relative to expectations for conventional quarkonia. 

In addition to this work on radiative transitions involving hybrids, another very interesting new development pursued 
by the JLAB LQCD group is the determination of the composition of the various states observed on the lattice, through 
the couplings of these states to a range of sources (such as radially and orbitally excited QQ operators, as well as 
hybrid operators) [l^. Of course in physical resonances the various basis states mix whenever allowed, and the 
mixing might a priori be large. However, in their test cases of higher-mass mesons (with a moderately large quark 
mass, between s and c) the JLAB group finds evidence that the basis states are not especially strongly mixed, so 
that dominantly radial, orbital and hybrid excited states can clearly be distinguished, and the pattern of states is 
broadly consistent with naive quark model expectations for the radial and orbital levels. This multisource approach 
to understanding resonances in LQCD is clearly very promising as a technique for improving our understanding of 
the nature of resonances in QCD. 

III. LOOPS IN CHARMONIUM: UNQUENCHING THE QUARK MODEL 

One of the long-standing mysteries in QCD is why the naive valence quark model describes the observed hadrons 
as well as it does. A simple quark potential model, even a nonrelativistic one, with linear scalar confinement, OGE 
forces, and a "constituent quark mass" describes the spectrum of both mesons (as qq states) and baryons (as qqq 
states) rather well, and also gives reasonably good predictions for EM, weak and strong transitions between these 
states. (For cc see for example Refs. [17r.lft |.) 

The fact that this type of model works reasonably well in describing hadrons is notoriously to justify. One might have 
expected various corrections to this valence quark model, such as relativistic quark motion, the effects of the gluonic 
degrees of freedom, and mixing between basis states (including higher Fock basis states) to completely invalidate this 
model. It is well established nonetheless that this model gives a reasonable first approximation to the spectrum of 
hadrons, and in heavy quark systems may even now be the most accurate theoretical approach available. This is 
one of the principal reasons for the great interest in the light charm-strange mesons D*q(2317) and Dsi(2460); they 
provide an example of a case in which the valence quark potential model clearly fails to predict masses accurately. 

Heavy quark hadrons provide attractive systems for the study of the reasons for the success of the naive nonrela- 
tivistic valence quark model, since one might expect these corrections to the quark model, such as relativistic effects 
and configuration mixing, to be less important. One specific aspect of configuration mixing that has seen an increase 
in recent interest is mixing of higher Fock space states, in particular the effect of an additional light qq pair in heavy 
meson systems. In one class of model it is assumed that the important mixing is between the quark model QQ basis 
state and the continuum of two- meson open flavor color singlet basis states, QQ o {Qq){Qq). This involves phe- 
nomenological but well-studied strong decay amplitudes, such as are usually described using the ^Pq model. This type 
of interaction, when treated to second order in the decay amplitude, predicts the mass shifts and level of (QQ) ^ (QQ) 
configuration mixing due to loops of open-flavor mesons, such as J/ip — )■ DD — >• J/i/j and ^Si {QQ) Di {QQ). 

Unfortunately, when numerical studies are carried out with physically realistic QQ -H> {Qq){Qq) channel couplings, 
one finds that the mass shifts due to specific individual {Qq) {Qq) channels are so large (and vary considerably from one 
QQ state to another that it is difficult to see how the naive valence quark model, which neglects these channel-coupling 
effects, could ever have been successful. 
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This conundrum may have been resolved recently in a study of coupled-channel effects in charmonium reported by 
Barnes and Swanson ^20|], who used the ■^Pq strong decay model to describe the valence-continuum couplings. This 
reference found that although individual intermediate open-flavor continuum channels made very different contribu- 
tions to the mass shifts of individual valence states, when a loop sum was carried out over complete spin multiplets 
(all j meson states from an n^i (Qq) + h.c. multiplet), the total mass shift was found to be quite similar for every 
QQ state within an N,L multiplet. This is evident in Table Hill which shows the DD, DD*, D*D*, (cs) analog, 
and summed mass shifts for the IS, IP and 2S cc states. Clearly the individual state mass shifts due to any specific 
intermediate channel are strongly state- and channel-dependent, but the summed mass shifts are rather similar, and 
could have been largely absorbed in a redefinition of parameters. 



Bare cc State 




Mass Shifts by Channel, AMi (MeV) 






Multiplet State 


DD 


DD* 


D*D* 


D,D, 


d,d: 


d:d: 


Total 


p - 

-1- cc 


IS j/iPii-'Si) 


-23 


-83 


-132 


-21 


-76 


-123 


-457 


0.69 


r?c(l'So) 





-114 


-105 





-106 


-98 


-423 


0.73 


2S V'(2''Si) 


-27 


-84 


-126 


-19 


-70 


-113 


-440 


0.51 





-118 


-103 





-102 


-94 


-544 


0.61 


IP X2(1'P2) 


-40 


-105 


-144 


-33 


-88 


-111 


-521 


0.49 


Xi(l'Pi) 





-127 


-148 





-90 


-130 


-496 


0.52 


Xo(l'Po) 

/lc(l'Pl) 


-57 





-196 


-34 





-172 


-459 


0.58 





-149 


-130 





-118 


-107 


-504 


0.52 



TABLE III: Mass shifts and continuum mixing induced in charmonium states by open-charm {cq){qc) meson loops, as found 
by Ref.[23|. Note that the mass shifts due to individual vary widely, but the sum is roughly state independent. Also note 
(final column) that the resulting physical states have rather large continuum components. (Pec = |{'I'|cc)p is the overlap of the 
physical state with the original cc valence state.) 

Even more remarkably, given certain rather mild constraints on the members of these multiplets (identical spatial 
wavefunctions within each N,L QQ multiplet, and separately within each n, £ Qq multiplet, initially degenerate 
valence states within each multiplet, and '^Pq decay couplings), Ref. (2Qj found that this was an exact result: An 
initially degenerate valence multiplet remains degenerate when coupled-channel (hadron loop) effects are included. 
(The mass shift is the same for all QQ states in a multiplet.) This explains how loop effects on the hadron spectrum 
could be both large and hidden; they were subsumed in an overall mass shift that is approximately state independent. 

One can also prove two other related results, regarding loop- induced mixing of different valence basis states through 
loops (which cancels in this limit) and strong widths (which are all identical in this limit, which imposes identical 
phase space for decays). This leads us to the "Three Lav^rs of Loopotics" (with apologies to Asimov [21|), which are 
that given the constraints specified above, the following results hold for the QQ states {A} in a given N,L multiplet: 

1. The mass shifts due to hadron loops for all states {A} in a given N,L multiplet are equal. 

2. Their strong (open-flavor) total widths are also equal. 

3. The configuration mixing amplitude a/i between any two valence basis states i and / due to 
hadron loops vanishes ii Li ^ L f or Si ^ Sf. 



These results have since been extended to more general decay models by Close and Thomas 1231 . We also note 
that a similar result regarding common loop mass shifts was previously reported by Tornqvist [23|. Recent studies 
of loop effects on charmed hadrons or charmonia include work by Baru et al. |24l |. Gamermann et al. 25|, Guo, 
Krewald and Meissner [2^, Huang and Kim [23|, Simonov and Tjon [2^, and van Beveren and Rupp f2£|. Since the 
continuum components of physical hadrons are typically found to be relatively large, it would be interesting to identify 
a "smoking gun" experimental measurement that is sensitive to loop components, such as a radiative transition in 
which the photon couples much more strongly to a Qq loop meson than to the heavy-quark QQ valence state. 
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IV. PANDA CROSS SECTIONS AND BES 

The PANDA project at GSI [sO] proposes to produce charmonia, and ultimately J^'-^-exotic charmonium hybrids, 
using pp annihilation. Although pp annihilation has been used previously to study charmonium, by E760 and the 
follow-on experiment E835 at Fermilab, these earlier studies used direct s-channel annihilation [pp cc ^ hadrons), 
which only produces non-exotic J^'-' states. Exotics will instead require associated production, in which a J^'-^-exotic 
charmonium hybrid iJcc (for example) recoils against another hadronic system, such as a meson m; pp ^ m + H^c- 

The PANDA project will only be feasible if these cross sections are sufficiently large. It is therefore somewhat 
unsettling that very little is known about these cross sections experimentally; the entire (published) data set on these 
associated charmonium cross sections consists of a few cross section measurements for the process pp — >■ ir'^J/ip near 
the he mass, which were taken by E760 and E835 for background estimates. (Data from E760 and E835 are discussed 
in Ref.jsij, which however does not quote acceptance-corrected cross sections for this reaction.) These earlier Fermilab 
measurements do provide a cross section scale (~ 100 pb), although we have no indication from this data of how these 
associated charmonium cross sections might vary with the choice of the light meson m, the charmonium state, or the 
i/s of the reaction. 

There have been several calculations of these cross sections using results from related processes, notably charmonium 
decays of the type (cc) — ppm, which have now been observed in several channels, and are related to the PANDA 
reactions by crossing. The earliest of these calculations (Gaillard et al. [s^, from the LEAR era) introduced a 
hadronic pole model, and discussed angular distributions in pp — >■ ir'^J/ijj at one energy. Lundborg et al. [s^ next 
derived a crossing relation, and used T{J/tp — >• ppn^) to predict a{pp — > ir'^J/ip) in a simple constant-amplitude model; 
the results were about a factor of 2-3 larger than the E760 cross section. The Gaillard et al. hadronic pole model for 
pp — > n'^J/ip was next generalized to other channels by Barnes and Li [34] . who reported cross section results for a 
wide range of PANDA {pp m + {cc)) reactions. Other interesting related results include evidence for a Pauli term 
in the J/ippp hadronic vertex [35*1 (suggested by BES data for the angular distribution of e+e" — >■ J ftp — >■ pp) and 
(note added in proof) the possibility of extracting NNm meson-nucleon strong couplings directly from the decays 
(cc) -i'ppm at BES [S^. 
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FIG. 3: Theoretical {pp — )■ m + (cc)) and experimental {pp — )■ vr" J/i/), from E760 and E835) PANDA type cross sections, from 
Ref.m. 

A comparison between theory (Ref.flHl, for a range of pp — >■ m -I- (cc) cross sections), and experiment (E760 and 
E835, for pp ■n^J/tp only) is shown in Figl^l Note that the cross sections for producing the 77c and xo are much 
larger than the other cases considered, due to their large {cc)pp couplings, which are inferred from the (cc) — ^ pp 
partial widths. This may imply that cc states couple most strongly to pp through gg intermediate states, so that 
PANDA would find the largest associated production cross sections for charmonia or charmonium hybrids having gg 
quantum numbers. 

V. DEEP CHARM MOLECULES 

Another topic that has been affected by the renaissance in charm spectroscopy is the subject of hadronic molecules. 
The discovery of the X(3872) in J/V'tt+tt" in B decays by Belle [37j has given us our strongest candidate for a 
meson-meson molecular state, since its mass and quantum numbers are consistent with a weakly bound D^D*^ {h.c. 
implicit) pair, but are inconsistent with expectations for a conventional cc state. Theoretical models suggest that this 
J^^ = 1++ D^D*^ system experiences an attractive force due to one pion exchange that is just strong enough to 
form a weakly bound state (perhaps also requiring additional, subdominant interactions) . 
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Although the system is an S-wave hadron pair, which is clearly favored for binding, Close and Downum [38| recently 
made the interesting observation that the two hadronic vertices involved in this one pion exchange require a P-wave 
Dtt pair. (e.g. DD*{\-) ^ D[D{Q-)'k{Q-)]p [D{Q^)tt{Q-)]pD D*{1~)D.) They note that processes in 
which the pion is emitted and absorbed in an S-wave meson-pion system should experience much stronger forces 
in general, and (when attractive) may form more deeply bound charm meson molecules. This S-wave requirement 
suggests consideration of molecules of meson pairs with the same J and opposite P, such as D*Di. This idea has been 
extended to other heavy flavor systems by Close, Downum and Thomas [39|. 

The binding energies calculated for the S-wave pion systems are strongly dependent on the short-distance behavior 
assumed for pion exchange, so precise values for the predicted binding energies cannot yet be quoted; Close and 
Downum estimate that these "deep charm molecules" might be ca. 100 MeV below threshold, hence near 4.2 GeV for 
D*Di, which suggests that the Y(4260) be considered a candidate. A search for DDtt-k-k decay modes of the Y(4260) 
is motivated by this possibility, since a D*Di molecule would likely decay significantly through constituent decay, 
D* — >■ Dtt and Di Dtttt. (Again an antiparticle h.c. label is implicit.) 

VI. CHARMISCELLENY 

We next consider two short items in charmonium spectroscopy which do not fit into the major categories discussed 
above, relating to the ■(/'(4415) and the Z(3930). 

The ^/;(4415) was for decades the highest mass charmonium state known, and in cc potential models is most often 
given the assignment 4'^Si. Although this state appears as a clear peak in i?, until very recently no exclusive decay 
modes had been identified, so the PDG (also for decades) listed this state as decaying dominantly to "hadrons", 
which one might have considered a safe assumption. This situation has finally improved with the observation of 
an exclusive -0(4415) decay mode by Belle fio']; it is seen in (2460)1?, with a combined branching fraction for 
5(0(4415) Dl?^(2460)) x B(D^{2460) D-tt+) of w 10 - 20% (depending on the 0(4415) parameters assumed). 
On considering other charge states and additional channels such as 1)5(2460) — >■ D*tt, this appears roughly consistent 
with the theoretical prediction of Ref.^lB], which anticipated that this was the second largest mode of a 4'^Si cc 
■0(4415), with a branching fraction of ~ 30%. As the largest mode is predicted to be DDiinarrow) (in a relative 
D-wave!), with a sa 40% branching fraction, a search for this mode is clearly of interest. 

And finally; the Z(3930) seen by BABAR [41] in two-photon fusion, 77 — > Z(3930) — > DD^ is now widely accepted 
as a 2^P2 cc state "Xc2(3930)", due largely to the agreement with the expected angular distribution for J=2, and to 
the plausible mass for a 2P cc state. It would be edifying to complete this picture through the identification of this 
state in 77 Xc2(3930) — >■ DD*; DD* is the only other open-charm mode available to this state, and the predicted 
relative branching fraction to DD* [19] is a reasonably large i?(xc2(3930) DD* /DD) « 0.35. 

VII. CONCLUSIONS 

In this invited talk I have reviewed some of the more recent developments in the theory of charmonium and related 
states; specific topics discussed included double charmonium production, LQCD and charmonium, open-charm hadron 
loops, charmonium production cross sections at PANDA, deeply bound charm molecules, and two additional (recent 
and possible future) measurements relating to cc candidates. It is evident that the spectroscopy of charmonium and 
related states is a remarkably active and interesting field of research in hadron physics. With the development of the 
BES-III and PANDA facilities, this happy situation should continue for many years to come. 
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